The second critical point for the Perfect Bose gas in quasi-one-dimensional traps 
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We present a new model of quasi-one-dimensional trap with some unknown physical predictions 
about a second transition, including about a change in fractions of condensed coherence lengths due 
to the existence of a second critical temperature T m <T C . If this physical model is acceptable, we 
want to challenge experimental physicists in this regard. 
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1. It is well known since 1925 jl[ for an ideal gas of iden- 
tical bosons of mass m in a three dimensional cubic box 
A = L 3 described the grand canonical ensemble (V, T, p) 
with T the temperature and p the chemical potential, 
there exists a critical density of thermal gas saturation 
of p c = £(3/2)/A| where Ap = h^/2-K(3/m is the thermal 
length and (3 = 1/fcsT . The hypothesis of A. Einstein 
back in 1938 [H] by F. London assumed that for a par- 
ticle density p above the critical density p c surplux the 
particles p — p c starts in the ground state mode k = . 
Thus we note that p k =o = p — Pc for P > Pc- But this 
argument macroscopic population of the ground state de- 
pends drastically on the geometry of the box. This is only 
first discovered H. Casimir Q in 1968 for a certain geom- 
etry anisotropic and what has motivated the work of M. 
van den Berg, J. Lewis and J. Pule from 1981 [H Q- 
These are generalized the concept of Bose-Einstein con- 
densate defining density in thermodynamic limit for an 
ideal Bose gas in a rectangular box any A = L x x L y x L z : 



Pa = lim 



lim 



\\k\\<v 



Pk = P 



(1) 



This definition includes all possibilities from a conden- 
sate, given here from Q and revisited in reference to 
new terminologies 

(i) The condensate conventional or usual [|| , [§| formed 
by the ground state (type I 0]). 

(ii) The fragmented condensate [U, @ if the general- 
ized condensate is distributed over a set (finite / infi- 
nite) modes macroscopically occupied in an energy band 
near the fundamental mode ie No = ^2 k<k Nk, with 
N k = O(N) (Type I / II @) 

(hi) The quasi-condensate [8J, |9j if the condensate 
is distributed over a generalized set of modes occupied 
mcsoscopic in an energy band near the fundamental mode 
ie = E fc<fc Nk, with N k = 0(N S ), S < 1 (type III 

In 1983 [5J, van den Berg proposed a box model three- 
dimensional anisotropic exponentially in both directions 
to define a second critical density p m separating a regime 
of condensate generalized type III (for p c < p < p m ) 
towards a regime of generalized type I condensate (for 
P > pin )■ Lately, we have reexamined this model 



flfjj questions that we call two-dimensional and we have 
shown that this a transition between two regimes: quasi- 
condensate and coexistence between a regime of con- 
ventional condensate and quasi-condensate, notably by 
showing that the second critical density p m corresponds 
to a saturation density of gas near condensation similarly 
to the saturation of thermal gas for p c . In addition we 
have calculated the second critical temperature and de- 
termined the fractions similar condensed changed. Then 
we calculated the effects of this transition on the coher- 
ence length for these two regimes of condensate general- 
ized. 

However this transition in the model boxes van den 
Berg is not the case for almost one-dimensional (quasi- 
1D). Our finding presented here is the existence of this 
transition between two regimes for the ideal gas of bosons 
in a harmonic trap quasi-lD. The purpose of our letter 
is to show where is the transition between a quasi-BEC 
and BEC conventional find the geometric model for the 
trap to his observation and to characterize this transition 
(Fractions condensed, coherence lengths). This theoreti- 
cal prediction is interesting because we have seen in the 
experimental 11 1 and theoretical articles [lij], 13 1 the 



existence of a second critical temperature called for 
quasi- ID environments, but calculated and explained dif- 
ferently from ours. 

2. Consider an ideal gas of bosons in a harmonic trap 
characterized by three-dimensional pulsation ui x 
where external potential is given by: 
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1 



1 



V ex t{r) = -mw x x + -mu; y + -mu z z 



The Hamiltonian for a particle is given by [14j : 
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= - — A + V ext (r) 
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The energy levels arc then given by: 



111 

e S:c s y s z = fr^z{s z + + ^z(Sz + jj) + tWz{Sz + g)- ( 4 ) 

Let p, the chemical potential of this system in the 
grand canonical ensemble, its value must be strictly lower 
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than the ground state to ensure the convergence dis- 
tribution the Bose-Einstein statistics showing the num- 
ber average boson in each of the modes of expression 
N SmSy s z = (exp s* ~ A 1 ) - Thus we can de- 

fine the effective chemical potential relative to ground 
level as /x = eooo — A/i. When the potential A/x > 
tends to 0, the thermal gas is saturated and we obtain 
a critical number of particles when ljq = (u^u^u^) 1 / 3 
tends to 0: 
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The nature of the spectrum for this trap is different 
from that of a particle in a box, that is why the type of 
geometry for the second transition will be different. If 
we consider that this trap is anisotropic with a quasi- ID 
type "cigar" in the axis z, we must have the condition 
lu z -C with lu x = Lo y = ui±. Usually 12 1, 13 1 is 

taken as a model lj±_ = \lj z with the condition A > 1. 
Thus we have an anisotropy level of modes in each direc- 
tion for hio z <C hu)±. If in this case we have an effective 
chemical potential hw z <C A/i -C Swi then we can cal- 
culate the number of atoms in a condensate generalized 
(quasi-condensatc) in an energy band near the fundamen- 
tal on the axis z: 



N qbec = N > 

s=(0,0,Sg) 
1 
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hf3uj z 



bx\J3An] =N-N C . (6) 



Thus for a macroscopic number of particles in 
the condensate, we must have condition A/i = 
p-i e -K0w z {N-N c ) But tllere ig anot h er condition 

for completing the energetic part along the axis z because 
the chemical potential A/i determines the size of the en- 
ergy band, decreases when the number of particles N 
increases and therefore a second critical number of par- 
ticles noted N m , the condition Huj z <C A/i is false. For 
this model unlike gas trapped in a box, there is no global 
notion of " density" of particles as the containment takes 
place throughout the three-dimensional space R 3 . How- 
ever we can introduce a scaling for the number of parti- 
cles n = Nujq which converges when N tends to infinity 



and when ujq tends to 
hw ~ A/i is a f3hu> z ~ e~ 



14] 



We then find that when 
-n c )/u>l _ Tq resolye thig 



condition, we introduce the model of anisotropic quasi- 

2/2 

ID trap following oj z = wie^^'" 1 where lo c > is a 
pulse of anisotropy exponentially. Thus there is a limit 
to the number of particles in the condensate distributed 
over the energy axis z: 



N m = N C + 



3 ' 



(7) 




Temperature 



FIG. 1: Illustration of phase diagram, see (0, J7]l 



Now if we go beyond this critical number N„ 
complete the ground state: 



we shall 



N bec = N S=Q = 



1 



1 



= N-N„ 



(8) 



and thus the chemical potential is A/i = l/(f3(N — N„ 
and the value of the quasi-condensatc is saturated: 



)) 



1 



Qbec 



npcj z 



ln[/3A/i] = N rn - N c 



since N — N c = (n — n c )/ojQ and thus uj^Nq^c ~ 

Note that this model of exponential anisotropy is nec- 
essary to limit "thermodynamics", ie when uj± tends to 
0, to define the second critical number of particles N m 
which is the same order as the first critical number N c 
(proportional to uJq 3 ). For finite pulses in a given exper- 
imental situation, wc must calculate the parameter oj c 

2/2 

depending on the aspect ratio A = e UJ ^' UJ± and see if the 
ratio u) c /uj± is sufficiently large to consider that one is 
in an exponential phase. Usually if one considers a sys- 
tem where very anisotropic A 3> 1 to consider this point 
of view that the system is exponentially anisotropic, we 
have to have ln(A) ^> 1 . 

3. In experiments with BEC, it is important to know 
the critical temperatures associated with corresponding 
critical densities. The first critical temperatures: T c , for 
a given particles number iVis given by solutions of this 
equation: N = N c = C^V^^o) 3 - Then wc get: 



haj 



fcsC(3)V' 



-TV 1 / 3 



(9) 



By the expression for the second critical particles num- 
ber one gets the following relation between the first and 
the second critical temperatures because T m are given by 
solution of this equation N = N m = N c + uj^/(h(3ujQ): 



T 



(10) 
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4. Another physical observable to characterize this sec- 
ond critical temperature is the correlation function con- 
densate coherence length. In quasi-lD harmonic trap, 
the correlation function is locally defined by: 



(16) 



FIG. 2: Curve illustrating (|lip as a function of Jr-. 



where the wave function (f> s (r) are the well known three 
dimensional hcrmitc functions: 




FIG. 3: Illustration of condensate fractions (|13p ( — ), (|14[) (- 

- -), m (- - -) with a = § 



Here r = Wcft/fc^^) 1 / 2 is "effective" temperatures re- 
lated to the corresponding geometrical shape. 

We can compute the relative difference between both 
critical temperature due to (fT0|) , see figure [2] 
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(11) 



defined by = ^ and with: 



t-2/3 



(12) 



where Q(x) = |(1 + Jl + ^£ 2 )- 

The temperature dependence of the "quasi- 
condensate" NQb ec is, see figured] 



N, 



Qbec 



N 



T rp 
J^ 1 



T m <T < T c 
T < T m . 



(13) 



The corresponding ground state conventional BEC is, see 
figure [U 







T < T < T 
± m — ± — ± c 

T < T 



N -|1-(X)3^ 
and for the two coexisting condensates one gets: 



(14) 



N — N c N N Qbec + N bec 
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(j) s (r) = (j) Sx {x)<j) Sy {y)(j) Sz {z), 
with for v = x,y,z (r x = x, r y = y, r z = z): 

1 



4>s v ( r u) 



-.H Sv (r v /L v )e 



with c Su = l/v s u \2 s » y/n, v = x,y,z and where the 
effective confinement lengths L v , v — x,y, z are given by: 

L v = J Notice that the limiting diagonal function 
n(r) = a(r, r) is local space particle density at point r. 

To compute the correlations in the direction z, in- 
duced by the quasi-condensatc, we have to consider 
r = r-Ar/2 r' = r + Ar/2, where Ar/2 = (0,0, Az) and 
integrate r over the space to keep the effect of condensate 
through the trap, it give the mean value of correlation for 
the distance Ar: 

(g)(Ar) = / d 3 rg(r- Ar/2, r + Ar/2). 
= (g Nhec ) (Ar) + (g Qbec ) (Ar) + (g bec ) (Ar) (17) 

where correlation function is divided in three parts, 
each part is the truncation of the sum ([TB")) : the non- 
condensate correlations (gNhec) = J2 S >o s >o s > * nc 
quasi-condensate correlations (gQbec) = =0 s =0 s 
and the condensate correlations (g bec ) = (correlation 
of the ground state). By formulas (f!6| and (fTT]) and 
by mathematical assumptions [l6|, |l7l | for the non- 
condensate part, we get, see figure |U 



(gNhec) (Ar) ~ / dx dy dzS~] 
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(18) 



which decrease exponentially fast for /i < 0. The quasi- 
condensate part, for T < T c , by the scaling argument of 
localization of the scale of energetic momentum of quasi- 
condensate given by Afi, and by mathematical assump- 
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between these two critical points. The next work will deal 



FIG. 4: (a) Allure of correlation functions (|19p with oj = 1, 
pti 2 /2m = 1 and /3A^ = 1, 0.1, 0.01, 0.001, 0.00001, 0.0000001 

respectively curves ( ), ( ), (...) and ( — ). Notice that 

(...) and ( — ) are combined, what can be interpreted as a 
saturation of the coherence of the thermal gas. (b) Allure 
of correlation functions (|18[1 with ujo = 1, ph 2 /2m = 1 for 
/3A/x = 0.001,0.0001,0.000001 respectively ( ), ( ), (— 




curve illustrating the exponent (|21[) as a function 
- |. This gives us a relative difference 
between the two critical densities of 0.587, see relation [Til and 
figure [2] 



FIG. 5: 

of Sr and with — = 



tions [16| , [17j , is given by, see figure ID 



(g Qbec ){Ar) ~ / (izVi 

1 3 = 1 



x exp[— jm/3(u) 2 z) /3A/x] exp[— 2m 
= E^7T 1 -exp[-2 m - (Az)2A ^ 



(Ar) 2 /3A/x 



(19) 



The condensate part, for T < T m is given by: 

(6W>(Ar) = / dr4> 000 (Ar)(j) 000 (-Ar)N be 

JR 3 



Az 2 

= N bec exp (-pj-) ^ N bec 



(20) 



We can defined the coherence length of the quasi- 
condcnsate as the maximal scale of length of Az such 
that the correlations are non- negligible when uj z , uj±_ 
tends to zero. Then by a scaling argument, see equa- 
tion (|19[) we obtain L c = y h/mu>±(3AiJ,, where A/i is 
give n by p- i e W<»z{N-N e ) for N( _ < N < N ^ Thus 

L C ^J mujj_/h = (lu^/uJz) 1 ^, where the exponent j(T) = 
hPujl{N ~ N c )/lu 2 for N c < N < N m and 7 (T) = 
h0wf(N m - N c )/uj 2 for N > N m . Using relations © 
and (fTTJ|) between N, N c , N m and T,T C , T m , we can find 
the temperature dependence of the exponent: 




1, T < T„ 



(21) 



7. We propose an experimental investigation of a weakly 
interacting Bose gas in quasi-one-dimcnsional-harmonic- 
trap, with the opposite of Thomas- Fermi regime to verify 
experimentally the effect of geometry on the kinetic en- 
ergy and on the coherence properties of the condensate 



with the implementation of the interactions between par- 
ticles and its effect on the second critical t emp erature as 
it is done for the first critical temperature [15| . An other 
question for quasi-two-dimensional boxes is the validity 
of Hohenberg theorem 0] which forbidden the existence 
of the condensate on the ground state for homogeneous 
system. 

We arc thankful to Markus Holtzmann for useful dis- 
cussions about the exponential anisotropy effect. 
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